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Abstract 

Explicit codes are constructed that achieve the diversity-muhiplexing gain 
tradeoff of the cooperative-relay channel under the dynamic decode-and- 
forward protocol for any network size and for all numbers of transmit and 
receive antennas at the relays. 

A particularly simple code construction that makes use of the Alamouti 
code as a basic building block is provided for the single relay case. 

Along the way, we prove that space-time codes previously constructed 
in the literature for the block-fading and parallel channels are approximately 
universal, i.e., they achieve the DMT for any fading distribution. It is shown 
how approximate universality of these codes leads to the first DMT-optimum 
code construction for the general, MIMO-OFDM channel. 

1 Introduction 

Cooperative relay communication is a promising means of wireless communication 
in which cooperation is used to create a virtual transmit array between the source 
and the destination, thereby providing the much-needed diversity to combat the 
fading channel. 

Consider a communication system in which there are a total of + 1 nodes 
that cooperate in the communication between source node S and destination node 
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D. The remaining {N — 1) nodes thus act as relays. We follow the literature in 
making the assumptions listed below concerning the channel. Our description is in 
terms of the equivalent complex-baseband, discrete-time channel. 

• All nodes have a single transmit and single receive antenna and are assumed 
to transmit synchronously. 

• The number of channel uses T over which communication takes place is 
short enough to invoke the quasi-static assumption, i.e., the channel fading 
coefficients are fixed for the duration of the communication, 

• We assume half-duplex operation at each node, i.e., at any given instant a 
node can either transmit or receive, but not do both. 

• The noise vector at the receivers is assumed to be comprised of i.i.d., circu- 
larly symmetric complex gaussian C7V(0, cj^) random variables. 

2 The DDF Protocol 

Under the DDF protocol, the source transmits for a total time duration of BT 
channel uses. This collection of BT channel uses is partitioned into B blocks with 
each block composed of T channel uses. Communication is slotted in the sense 
that each relay is constrained to commence transmission only at block boundaries. 
A relay will begin transmitting after listening for a time duration equal to b blocks 
only if the channel "seen" by the relay is good enough to enable it to decode the 
signal from the source with negligible error probability. We explain in more detail. 

2.1 Notation and Expressions for the Received Signal 

Initially, we will assume that each relay node has a single transmit antenna. The 
extension to arbitrary number of antennas is straightforward. We will similarly 
make the initial assumption that the destination node has a single receive antenna. 

It will be convenient at times to regai^d the source as the first relay, i.e., S = Ri 
and the destination as the (A^ + l)th relay, i.e., D = Rn+i- The notation below 
is with respect to a fixed channel realization that lasts for the B-b\ock duration. 

Let Xij{n), 1 < b < B, n = 1,2, ■ ■ ■ ,N denote the T-tuple transmitted by the 
nth node during the 6th block. Since all nodes do not transmit in all blocks, we will 
make the assignment a;^(n) = ip, where we regard if as the "empty" vector to han- 
dle the case of no transmission. In particular, the vectors b = 1,2, - ■ ■ , B 
denote the B successive transmissions by the source. 



2 



Let us assume that up until the end of the (6 — l)th block, we know which 
relays began transmitting and when. We will assume that once a relay has begun 
transmitting, it will keep on transmitting thereafter until the end of the Bth block. 
Let Xfc denote the set of indices of the relays that transmit during the kth block, 
k = 1,2, ■ ■ ■ , B. We will refer to 2f^ as the kth activation set. Clearly 

^1 = {1} 

Ik ^ Ik+u l<k<{b-2). 

We next proceed to determine for the relays not in Tt-i, whether or not the time is 
right for them to begin transmission during the 6th block. In other words, we will 
determine lb given {1k}^fZ}i- Since Ii is known, this procedure will allow us to 
recursively determine the activation sets Ik for all 1 < A; < i?. 

We will begin by first identifying the signal received by such a relay during the 
(fe — l)th block. Let Cb, 1 ^ t> < B, denote the size of i.e., 

\Ib\ = Cb- 

Clearly, 

1 = Ci < C2 < • • • < Cfe-i < iv- 

Let the elements ofl^, 1 < /c < (6 — 1), be given by 

Ik = {1 = mi, 7712, • • • , J. 
We use h{m,n) to denote the fading coefficient between the mth and nth nodes. 

[h{mi,n), h{m2,n), ■■■ , h{m(;^,n)] 
Xfc(mi) 
Xk{m2) 

_ Xk{mc_^) _ 

[y(k,i){n) y(k,2){n) ■■■ y{k,T){n)] 
[w{k,i)(.n) W(^k,2){n) ■■■ W(^k,T){n)] 

denote the received signal and noise vector at the nth node during the A;th block. 
Then we have 

l/kin) = hi{n)Xk + wi{n). 



Let n Ib~i and 

hkin) 

Xk 



Let 



llk^n) 
dkin) 
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Therefore the totality of the received signal at the nth node up until the end of the 
{b — l)th block is given by 



+ Wi{n) ■■■ 2«l-i(n)]. 



(1) 



Note that the vectors (n) as well as the matrices , 1 < / < 6 — 1 are in general, 
of different sizes. 



2.1.1 Signal at Destination 

Since D = Rn+i, by replacing n by {N + 1) and 6 — 1 by S in equation ([T]) above, 
we recover the expression for the received signal at the destination during the Bth 
block: 



\IA{N + 1) 



h%{N + l)] 



Xb 



(2) 



2. 1.2 Outage of Relay Node 

From ([T]), we note that the channel "seen" by the nth relay node over the course of 
the first 6—1 blocks is the MISO (multiple-input single output) channel character- 
ized by the matrix equation 

y = ilAin) ••• li-i{n)]x + w. (3) 

The nth relay node can only hope to decode reliably at the end of the (6 — l)th 
block if at that point, it has sufficient mutual information to recover the transmitted 
signal whose information content equals rB\og{p) bits. Here r denotes the mul- 
tiplexing gain, p the signal to noise ratio, and r\og{p) the rate of communication 
between source and destination [15]. If it does not have sufficient information, then 
we say that the relay is in outage. Thus the probability of outage Voax,n,b-i{f) of 
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the nth relay node at the end of the [b — l)th block is given by 

^out,n,6-l('^) 

= Pr (^(1 + I h^in) p) < r^^^^log(p)^ 

= Pr 1^(1 + p I biin) < log(p)^ . 

Under the DDF protocol, the nth relay node at the end of block b — 1 uses this 
expression to decide whether or not it is ready to decode. If it is ready to decode, 
then it will proceed to do so and then begin transmitting from block b onwards, i.e., 
n € Xft. 

2.2 Performance under the DDF Protocol 

In our analysis of the DDF protocol, we will make the assumption that if a relay 
does decode erroneously, then this error will propagate and cause the receiver to 
decode incorrectly as well. Thus the receiver at the destination will decode cor- 
rectly if and only if in addition to the receiver at the decoder, the receivers at all 
intermediate nodes that have participated in relaying of the transmitted signal have 
also decoded correctly. 

A lower bound on the probability of error of the DDF scheme can thus be 
derived by making the assumption that when the channel seen by a relay node 
is not in outage and the relay proceeds to decode the signal transmitted by the 
source, it will do so without error. Under this condition, the error probability of the 
DDF scheme, will be lower bounded by the probability of outage of the channel 
dl]), seen by the destination. In Section 13.21 we will construct codes whose eiTor 
performance at large SNR is equal to this lower bound, thereby establishing that 
this lower bound is indeed the error probability associated with the DMG tradeoff 
of the DDF protocol. 

Let 7 denote the vector composed of the (^^^) fading coefficients 

iu^ l<m<N, ] 

|Mm,n)| n>m, 2 < „ < (iV + 1), / 

ordered lexicographically. We will use T to denote the random vector of which 
7 is a realization. The activations sets Tk are clearly a function of the channel 
realization 7. Writing Ikil) in place of to emphasize this, let us define 

lij) = (Xi(7),--- Ml))- 
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Let A denote the collection of all possible activation sets. It follows that the error 
probability of the DDF scheme satisfies 



where 



7^(J) 
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T(7) =1 
^ + PT.ti\h^i{N + l) \A<r\og{p) 



2.3 Notation to Aid in Code Analysis 

Returning to the expression for the signal at the nth relay node up until the — l)th 
block in ([I]), we extend the vectors h^{n) and the matrices to be of equal size 
with a view towards the ST code construction to be presented in Section [l!2l 
The vectors 

{lkk{n) I 1 < A: < 6 - 1, 1 < n < (iV + 1)} 

will be extended by zero padding, while the matrices X]^, I < k < b — Iwill be 
padded with arbitrary row vectors. The extra row vectors can be chosen arbitrarily 
since the extended matrix X^ will be left multiplied by row vectors /i^ (n) having 
zeros in the locations corresponding to the indices of the row vectors where padding 
of the matrix Xj. takes place. 

We thus define, for 1 < A; < 6 — 1, 



[hkil,n) hk{2,n) • • • /ifc(iV, n)] 



where 



Also, let 



hk{m,n) 



Xh 



h{m, n) m £ 2k 
else. 



^fc(2) 



where 



Xkim) 



Xk{m) m £lk 
arbitrary n-length vector else. 
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In terms of the extended vector and extended matrix notation, the received 
signal at the nth relay node, n 1b~i and the destination can respectively be 
re-expressed in the form 



[k{n) 



Xb- 



w 



b-l 



in)], 



(4) 



Xi 



Xb 



+ [nAiN + l) ■■■ + 



(5) 



In this representation, all vectors hi{n) are of the same size, (1 x T). The same 
comment also applies to the matrices 1 < / < 6 — 1, which are of size (A^ x T) 

As will be shown in Section |4]below, ST codes that are approximately universal 
for an appropriate class of block-fading channels will be the building blocks of 
codes for the DDF protocol that attain the DMG performance of this channel. For 
this reason, a discussion on the block-fading channel is presented in the next two 
sections. 



3 The Block-Fading Channel 
3.1 Outage Probability 

Consider the block-fading MIMO channel with rit transmit and receive antennas 
and B blocks, characterized by 

y_f, = HbX^ + Wf,, l<b<B. (6) 
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Thus each matrix Hi, is of size (n^ xnt). The probability of outage of this channel 
is given by 

B 

Pr{y^logdet{In^+ pHbHl) < rBlog{p)) 



Pout(r) 



b=l 



= Pr(log det(/B„^ + pAhA^h) < log(p)) 
= Pr(log det(/Bni + /oA^Ah) < rB log(/>)) 
where p is the SNR and where Ah is the {Brir x Brit) block diagonal matrix 

Hi 

H2 



A 



H 



Hb 



In the above, 
example, y = 



and <, > corresponds to exponential equality and inequality. For 

X. Let q = Uf B and let 



is used to indicate that lim — 



Ai < A2 < 



log('i/) 
log( 

< A„ 



be an ordering of the q eigenvalues of A^^Ah. Note that if < nt, then 



Ai 



A, 



A 



{nt~nr)B 



0. 



(V) 



(8) 



Let 6 
by 

Then 



{[nt — nr]B)'^ where (x)+ denotes max{x, 0}, and let the ai be defined 

6+l<i<q. 



A,; 



i^out(r) 



P 



Pr( X] (1 ~ "»)^ < ^^)- 



=5+1 



We will now proceed to identify a ST code in the next section. Section [ 
that is approximately universal for the class of block-fading channels, i.e., a code 
that achieves the D-MG tradeoff of the channel model in (Q) for every statistical 
distribution of the fading coefficients {[Hh]ij}. 

Similar construction of codes for such a setting have previously been identified 
in [20,26] and independently in [21]. We adopt the code-construction technique of 
these papers for the most part, although the construction presented here is slightly 
more general, for example, we permit the individual block codes to be rectangular 
and offer flexibility with respect to number of conjugate blocks employed. Most 
importantly though, our proof will establish the new result that these codes are 
approximately universal for the block-fading channel and parallel channels. 
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3.2 Approximately-Universal Codes for the Block-Fading Channel 



3.2.1 Constructing the Appropriate Cyclic Division Algebra 

Let T be an integer satisfying T > nt. Let m > Bhe the smallest integer such that 
the god of m, T equals 1, i.e., (m, T) = 1. Let K, M be cyclic Galois extensions 
of Q(^) of degrees m, T whose Galois groups ai^e generated respectively by the 
automorphisms (Ji, i.e., 

Gal(IEC/Q(z)) = <(t)i> 
Gal(M/Q(z)) = < fji > . 

Let L be the composite of K, M, see Fig{T] Then it is known that L/Q(z) is cyclic 
and that further, 

Gal(L/Q(^)) ^ Gal(]K/Q(i)) x Gal(M/Q(z)). 

Thus every element of Gal(L/Q(?)) can be associated with a pair o-J) belong- 
ing to Gal(]K/Q(i)) X Gal(M/Q(z)). Let (p, a be the automorphisms associated to 
the pairs ((/)i,id), (id, o"i) respectively. Then 0, a are the generators of the Galois 
groups Gal(L/M), Gal(L/]K) respectively. 




Figure 1 : Construction of the underlying cyclic-division algebra. 

Let 7 G K be a non-norm element of the extension L/IK, i.e., the smallest expo- 
nent e for which 7*^ is the norm of an element of L is T. Let z be an indeterminate 
satisfying z'^ = 7. Consider the T-dimensional vector space 

D = {z'^~Ht-i e z^-^£t-2 e • • • ^0 Ki g l}. 

We define multiplication on D by setting £iZ = za{ii) and extending in a natural 
fashion. This turns D into a cyclic division algebra (CDA) whose center is K and 
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having L as a maximal subfield. See [13, 16] for an exposition of the relevant 
background on division algebras. Every element x = z^~^£t-i + z^~'^£t-2 + 
■ ■ ■ + l^m D has the regular representation 



X 



£t-i cf{£t-i) 



(9) 



The determinant of such a matrix is known to lie in K. Given a matrix X with 
components Xi j € L, we define (t){X) to be the matrix over L whose {i,jy^ 
component is given by [^(X)]^ j = (l){[X]ij). Note that in this case, 



m—l 



m—1 



l[det{cP\X)) = ll^\det{X)) 



1=0 



i=0 
m—l 



n '^i(det(X)) 



i=0 



Hence if the elements ii underlying the matrix X are in addition, restricted to lie 
in the ring Cl of algebraic integers of L, then we have that 



m—l 



so that 



1=0 



m—l 



n det(<^^(X)) 



> 1. 



(10) 



3.2.2 Space-time Code Construction on the CDA 

Let X be the rectangular {iit x T) ST code comprised of the first nt rows of the 
regular- representations of the elements YlJ=o where ii are restricted to be of 
the form: 

T 

ii = ^^jjTji € >1qam 



10 



where {71, • • • , 7t} are a basis for h/K and where 



QAM 



{a + ib I \a\, \b\ < (M - 1), a,b odd} C 



denotes the QAM constellation of size AP. Note that as a result, we have ensured 
that ii € Ol- Thus each code matrix in A' is of the row-deleted form 



X 



£0 7(t(£t-i) 



7CJ 



T-li 



(11) 



Let S be the {Bnt x BT) ST code comprised of code matrices having the block 
diagonal form: 



X 



c^^~\X) 



, X 



where 6 accounts for SNR normalization. When this code matrix is in use, the 
received signal over the block-fading channel is given by 



B\ 



[Y1Y2 

This can also be expressed in the form 



[Wi W2 



Hb]S + 
• Wb]. 



(12) 
(13) 



Yi ' 




' Hi 




X 






Y2 








4>{X) 




W2 




= e 






+ 




_ 




_ 











,(14) 



in which the channel matrix is of block-diagonal form. This latter form is conve- 
nient when comparing the block-fading channel with the parallel channel. 



3.2.3 Proof of Optimality 

We will now show that the ST code S is approximately universal for the class of 
channel models described by From information rate considerations we must 
have (M^)™^' = p''-^^, i.e., = . Also, since O^M"^ = p, we have 

(7 — p mT . 
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We next examine the error performance of the code. Let 

H = [Hi H2 ■ ■ ■ Hb] 

and let AS" = Si — S2 where 5i , ^2 , are two distinct code matrices belonging to 
S. We have the following expression for the squared Euclidean distance: 



e^Tr{HASAS^H^) 



> 



ntB 
i=l 



by the mismatched eigenvalue bound (see [13,30]), where /xi > /X2 > • • • > ^J-mB 
are the ordered eigenvalues of AS'AS^. 

Let AX be the T x T matrix that corresponds to AX in the sense that it 
would have been the matrix obtained if in constructing the ST code X, we had not 
deleted the bottom {T — nt) rows from the regular representations of the elements 



Correspondingly, let AS be the {BT x BT) matrix 



AS 



AX 



lB-I 



{AX) 



By the inclusion principle of Hermitian matrices, (see Theorem 4.3.15 of [31]) 
the smallest Brit eigenvalues of AS'AS^ are term-by-term lai^ger than the corre- 
sponding smallest Brit eigenvalues of ASAS\ i.e., 

HBnt-k+l > VST-k+l, 1 < ^ < Bnt, 

where z^i > > • • • > t^bt ai'e the eigenvalues of AS'AS'^^. 

Next, let AS be the {mT x niT) matrix obtained by further extending A^ to 
include all m "conjugates" of AX, i.e.. 



AS 



AX 



Let z>j, 1 < i>i < mT, be the eigenvalues of AS'[AS']^ ordered such that 

= I'i, 1 < i < BT. 
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Then for every I < J < q = Brit, we have 
J 



where 



d'E > ^ ^q-ifJ-q-i 

J 



i=0 
J 



^ \-ii'BT- 



i=0 



J J 

1 r 1 



i=0 



^ pmT-(J+l) ^ ^■'^^ 



6j = mT{l-^) + J + l-mT-Y, 



J + l-rB 



mT 

J 

"q-i 



a, 



q-i 



i=0 



eta 

1=0 



i=q—J 

In this derivation we have made use of the fact that the product of the eigen- 
values is equal to the determinant and of the non-vanishing determinant property 
enunciated in (flOl ). We will now show that if the block-fading channel is not in 
outage, that for some J, 1 < J < (5j > 0. 

Suppose 

Q 

5^(l-ai)+ > (r + e)S. (15) 

i=l 
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Clearly for some i, ai < 1. Since ai > a2 > ■ ■ ■ > oiq, assume ai > a2 > 
• • • > Oiq-j-i > 1, and ai < 1, i > q — J. Then ( fT5] ) becomes 

1 

i=q~J 

and this causes the corresponding 6j > eB hence leading to negligible error prob- 
ability when not in outage. By taking the limit as e — > we see that the probability 
of error is negligible in the no-outage region. This proves that the space-time code 
X achieves the DMG tradeoff regai^dless of the statistical distribution of the fading 
coefficients that comprise the matrices {Hi^}, i.e., proves approximate universality 
of the constructed ST code. 

3.3 Analogous Results Hold for the Parallel Channel 

By parallel channel we will mean the channel given by 







' Hi 




" Wi ' 








S + 




Yb 




■•• Hb 




Wb 



in which the channel matrix is of block-diagonal form. Consider the {Bnt x T) 
space-time code 5pai- given by 





X 










< e 









which when used over the parallel channel leads to the equation below for the 
received signal at the receiver, 



Yi ' 




' Hi 




X 




' Wi ' 


Y2 


= 9 


H2 






+ 


W2 


Yb _ 




■■■ Hb _ 











Comparing this equation with the alternate expression for the block-fading channel 
given in (fT4l) we see that the expressions are identical. There is one important 
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difference though. In the case of the block-fading channel, a rate requirement of 
R bits per channel use translates into a space-time code S of size 2^^^ = p^^^ , 
whereas in the case of the parallel channel, the size of the corresponding ST code 
5par is required to be 2^^ = p^'^ . 

It follows from this that by replacing rB by r, one can similarly prove approx- 
imate universality of the code S^^v for the class of parallel channels. We omit the 
details. 

3.4 DMT-optimal Codes for the General MIMO-OFDM Channel 

The MIMO-OFDM channel can be regarded as a parallel channel in which each 
parallel block corresponds to a different subcaiTier and can thus be represented in 
the form: 

y. = OHiXi + Wi, 1<1<Q, 

where where Q is the number of OFDM tones or sub-carriers [25]. The matrices 
Hi are correlated in general, with a correlation derived from the time-dispersion of 
the original ISI channel. Since the code 5pai- is approximately universal, this means 
that the code Sp^- is DMG optimal when used over the MIMO fading channel. 
When the code 5par is used over the MIMO-OFDM channel, the received-signal 
equation will take on the form 







' Hi 




X 








= 6 






4>{X) 


+ 












_ (t)'^~^{X) _ 







DMT-optimal codes for the OFDM channel have previously been constructed 
in [20] and [22]. In [20], the authors provide a proof only for the case when the 
matrices Hn appearing along the diagonal are i.i.d. Rayleigh. The DMT-optimal 
construction in [22] is for the SIMO-OFDM case. We thus believe the results in 
this paper represent the first construction of DMT-optimal codes for the general 
OFDM-MIMO channel. 

4 Codes Attaining the DMG of the DDF Protocol 

We now show how ST codes constructed for the block-fading channel can be used 
to construct optimal codes under the DDF protocol. 
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We consider the DDF protocol as it applies to a communication system in 
wtiich there are a total of 1 nodes that cooperate in the communication between 
source node S and destination node D. 

As in Sections [T]2l under the DDF protocol, the source transmits for a total time 
duration of BT channel uses. This collection of BT channel uses is partitioned into 
B blocks with each block composed of T channel uses. Communication is slotted 
in the sense that each relay is constrained to commence transmission only at block 
boundaries. A relay will begin transmitting after listening for a time duration equal 
to b blocks only if the channel "seen" by the relay is good enough to enable it to 
decode the signal from the source with negligible error probability. 

Our coding strategy runs as follows. The role played by nt in the block-fading 
scenario is now played by the number which is the number of nodes in the 
network capable of transmitting to the destination. Let A! be the rectangular {N x 
T) ST code comprised of the first N rows of the regular representations of the 



elements ^ 



T-l 



where ^, are restricted to be of the form: 



T 

i=i 



Let V be the {BN x BT) ST code comprised of code matrices having the block 
diagonal form: 



V 



X 



<t>{X) 



X 



(16) 



where 9 accounts for SNR normalization. The code to be used then has the follow- 
ing simple description. The source 5 sends the first row of each of the matrices 
X, (piX), ■ ■ ■ , (l)^~^{X) in successive blocks. Let us assume that relay node Rn, 
2 < n < iV, is not in outage for the first time at the conclusion of the (6 — l)th 
block. Then Rn is ready to decode a the end of the b — 1th block. Thereafter, it 
proceeds to send in succession, the nth rows of the matrices (l)''{X), (l)''~^^{X), • • • , 
(j)B~i(^X). Thus the appended matrices Xi appearing in (01, con^espond to the 
matrices in (fT6l ). 

It is easy to show using the results stated earlier relating to the block-fading 
channel that this coding strategy ensures that whenever a relay node decodes, it 
does so with negligible probability of error. The destination eiTor probability is 
also similarly guaranteed to have error probability that is SNR-equivalent to the 
outage probability, thus proving DMG-optimality of the constructed ST code. 
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This follows since each relay node , n Ib~i "sees" a block-fading chan- 
nel (see (Ull and the coding strategy we have adopted ensures that the code matrix 
canying data from the nodes in T^-i is DMT optimal for the corresponding block- 
fading channel. A similar statement is true for the relay Rn+i that corresponds 
to the destination, since the corresponding channel equation is of the same block- 
fading form see 

4.1 Example 

We illustrate with an example. Consider a network in which there are a total of 
+ 1 = 5 nodes including source S = Ri and destination D = nodes. Let 
the block length T = 4 and the number of blocks B = 4. Let us assume that 
at the end of the first block, relay R^ is not in outage and therefore in a position 
to decode. Let us assume that relay R^ is ready to decode at the end of the 3rd 
block and that relay R2 is in outage throughout and thus does not participate in the 
communication. Thus 



here. In terms of the notation introduced in Section 12.11 the signal received by 
relays R2, R3, R^, R5 prior to decoding are given as follows. Decisions at each of 
the relays i^s, i?4, R^ ai"e based on the received signals 



1a 



{1}, 
{1,3}, 
{1,3}, 
{1,3,4} 



ri(3) 



/i(l,3)x*i(l) + 



bl(4)y*2(4)y*(4)] 



[/i(l,4) /i(l,4) /i(3,4) /i(l,4) /i(3,4)] 



x*(l) 

i£*2(l) 
^2(3) 

31^(3) . 



+ ki(4)«!^(4) 1^7^(4)] 
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[^1(5)^2(5)^3(5) yi(5)] = 

[/i(l,5) /i(l,5) h{3,5) /i(l,5) /i(3,5) /i(l,5) /i(3,5) /i(4,5)] 

r ^i(i) 1 

^*2(1) 

^l(l) 

xi(3) 

L 4(4) . 

+ [l^i(5) 2^2(5) li!*3(5) u;*(5)] 

respectively. 
Set 



X(l) : = 



" 3Ii(l) ' 




" 4 7C^(^3) 






^i(2) 




h ct(4) 






^i(3) 




^2 




7^3(^3) 


. a (4) . 




^3 aih) 


^'(^l) 





X(2) := </>{X(l)} 
X(3) := 0^{X(1)} 
X(4) := 

The corresponding extended vectors are given by 

[M(3)] = [Ml, 3) 0] 



^5(4) = [/i(l,4) 0] 
42(4) = [/i(l, 4) 0/1(3,4)0] 
^3(4) = [Ml, 4) 0/1(3,4)0] 



45(5) = [/i(l,5) 0] 

^2(5) = [/i(l,5) /i(3,5) 0] 

^3(5) = [/i(l,5) /i(3,5) 0] 

4*1(5) = [/i(l,3) /i(3,5) /i(4,5)] 
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and we have 



X{1) 



[yl(5)y*(5)y*3(5) y*(5)] 



[4*i(5)A2(5)43(5)Ai(5)] 



X{2) 



X{3) 



X{4) 



4.2 Extension to Multiple Antenna Case 

The extension to the case of multiple antennas at each relay node and at the desti- 
nation is straightforward. 

5 An Alamouti-based Code for the case of a Single Relay 

In this section, we provide a particularly simple code construction for the case 
when in addition to the source S and destination D, there is a single relay antenna 
i?2- The basic building block for the distributed space-time code is an Alamouti 
code. A separate proof is given here as the proof given in previous sections does not 
apply here, primarily because the Alamouti code is not an approximately universal 



5.1 Constructing the Appropriate Cyclic Division Algebra 

The cyclic division algebra is constructed along the same lines as before, with some 
differences, for example Q here plays the role of Q{i) earlier . 

Here the number of channel uses in each block equals 2, i.e., T = 2. Let 
m > B be the smallest integer such that (m, T) = 1, i.e., m is the smallest odd 
integer > B. SetM = Q{i) and let IK be a cyclic Galois extension of Q of degree 
m. Let the Galois groups of IC/Q and M/Q be generated respectively by the 
automorphisms (pi, ai, i.e.. 



Thus (Ji corresponds to the complex conjugation operator. Let L be the composite 
of K, M, see FigH 

Then it is known that 




Gal(K/Q) = 
Gal(M/Q) = < fji > . 



Gal(L/Q(j)) 



Gal(IC/Q(2)) X Gal(M/Q(i)). 



This result was first presented at a poster session in Allerton 2006 
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m 2 
<</>!> \ /<o-i> 



Figure 2: Construction of the underlying cyclic-division algebra. 



Thus every element of Gal(L/Q(^)) can be associated with a pair cr{) belong- 
ing to Gal(]K/Q(z)) X Gal(M/Q(z)). Let 0, cr be the automorphisms associated to 
the pairs (^i, id), (id, ai) respectively. Then 0, a are the generators of the Galois 
groups Gal(L/M), Gal(L/IC) respectively. Here again, a is the complex conjuga- 
tion operator operating on L. We note that a commutes with 0. 

Let 7 = —1. Then 7 is a non-norm element of the extension L/K, i.e., the 
smallest exponent e for which 7*^ is the norm of an element of L is 2. This fol- 
lows because the norm in L /K of any element in L is non-negative. Let z be an 
indeterminate satisfying = 7. Consider the 2-dimensional vector space 

D = {zii®lo I li G L}. 

We define multiplication on D by setting liZ = za{ii) and extending in a natural 
fashion. This turns D into a CDA whose center is IK and having L as a maximal 
subfield. Every element a = zii + £0 has the regular representation 



X 



4 -to 



(17) 



which we recognize as the familiar Alamouti code matrix. The determinant of such 
a matrix is clearly real and thus lies in K. Note also that the rows of every such 
matrix are orthogonal and hence the two eigenvalues of the matrix have the same 
magnitude. 

Given a matrix X with components Xij G L, we define <j){X) to be the matrix 
over L whose {i,jY^ component is given by [4>{X)]ij = 4>{[X]ij). Note that in 
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this case, 

m—1 m—1 



i=0 i=0 

e Q. 

Hence if the elements underlying the matrix X are in addition, restricted to lie 
in the ring Oj^ of algebraic integers of L, then we have that 

m—1 

Yl det{(j)\X)) e Z 

i=0 

SO that 

m—1 

I n det(</>^(X)) |2 > 1. (18) 

1=0 

5.2 Channel Models and Outage 

We assume as before that the total duration of communication is 2B channel uses, 
partitioned into B blocks of 2 channel uses each. Let us assume that at the end of 
the (6 — l)th block, the relay determines for the first time that it is not in outage 
and begins to transmit from block b onwards. The channel perceived by the relay 
antenna is given by 

y = hsrX + w 

where w is the usual additive noise. The probability of outage of this channel is 
given by 

Pr(log{l + p I ft„ |2)) < ^^'^ log(p) 
= Pr((l-A)-<^) 

where /3i is defined by 

The channel seen by the destination takes on the form 

y = \hsr ■ ^ • /isr hsr Kd ' J hgr Kd] x + W 
(b — 1) terms 2[B - (b — 1)] terms 
= h^X + W, 
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where we have defined 

h . — [^sr * * * ^sr ^sr ^rd ' ' ' hgr ^rd]- 



(19) 



5.3 DMT-Optimal Code Construction 

Optimal code construction proceeds as follows. The source transmits the B blocks 

0[Acl>{A) ■■■cl>^-\A)] 



in succession, where 



and where 



A = [io - 1,] 



where 7j is a basis for L/Q(z) and where lij € ^qam- 

The relay transmits from block h onwards and its transmissions are of the form 

e[(i)\c) /+i(C7) ••• <^^~^(C)] 

where 

C = [h t^]. 

The signal seen by the receiver at the destination is thus of the form 

A 



y = 9h 



+ w, 



with h given in (T9\ . 



22 



5.4 Proof of Optimality 

We will show as in Section|4]that when either the relay or the destination is not in 
outage, the error probability incuiTed by this code is negligible i.e., of order p^°° 
thus proving DMT optimality of the code. 

Note that from rate considerations, we must have 



(M 



2\2m 



P 



2tB 



5.4.1 Optimality in the Broadcast Phase 

Let 

^AQAM = {a + z6 I a, 5, even, < |a|, |&| < 2(M - 1)} . 

The Euclidean distance between the received matrices associated to code ma- 
trices Xi , X2 is given by 



> 



b-1 



^0) I' + I <P\-t^) n , ^^ G ^AQAM 



i=0 
b-1 

i=0 

rb-1 

4 = 



hsr 



1 

6-1 



1 



2 / p 



P 



(m-b+l) 
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Consider the no-outage region associated to rate r + e: 



b-1 



and it follows that the probability of error is neghgible for all e > 0. 
5.4.2 Code in the Cooperation Phase 

Here, by making use of the orthogonality of the rows of the Alamouti code, we can 
bound the minimum Euclidean distance as follows: 



4(AX) > 02|;^,,|2^0*(|4|2 + |_4|2) 

B-1 

+OH\hsd\^ + \hrd\')Yl 

i=b 

B-1 

i=0 i=b 
b-l 

= 0^ I Kd ? E '^'W + ? + I Kd ?] E '^'W- 



B-1 



B-1 



i=0 



Let us abbreviate and write 



hi 

h 



h2 = y/\ hsr p + I Kd 



[ hi - - J hi /i2 ■ - h2^ ] 

b— Iterms B— b+ltemis 



Then this squared Euclidean distance is also the squared Euclidean distance 
over the block-fading channel shown below: 



y = Oh 



kB-l 



(^) 



+ W. 



Note that since 



|2 
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both channels are in outage for precisely the same set of values of the fading coef- 
ficients hgd, hrd- On the other hand by our results in Section I3.2.3I the block di- 
agonal code appearing in the equation above has negligible eiTor probability when 
the channel is not in outage. This proves DMT optimality of the code in the coop- 
eration phase as well. 
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